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Abstract. - We study structural phase transition of polymer-grafted colloidal particles by Monte Carlo 
simulations on hard spherical particles. The interaction potential, which has a weak repulsive step outside 
the hard core, was validated with use of the self-consistent field calculations. With this potential, canon- 
ical Monte Carlo simulations have been carried out in two and three dimensions using the Metropolis 
algorithm. At low temperature and high density, we find that the particles start to self-assemble and finally 
align in strings. By analyzing the cluster size distribution and string length distribution, we construct a 
phase diagram and find that this string-like assembly is related to the percolation phenomena. The average 
string length diverges in the region where the melting transition line and the percolation transition line 
cross, which is similar to Ising spin systems where the percolation transition line and the order-disorder 
line meet on the critical point. 



Colloidal particles are solid particles made from polymers, metals, and other various materials. 
Latex is an example of such colloidal particles. They show characteristic phenomena, for example 
depletion interaction and glass transition. It is known that such colloidal systems can be modelled 
by a rigid spherical particle system [1]. The phase behavior of rigid particle systems shows a phase 
transition between the solid phase and the fluid one, which is called Alder transition [2, 3]. 

Since colloidal particles tend to aggregate and flocculate, polymers are often grafted on the col- 
loidal particles used in industry to prevent such flocculation. They are called polymer-grafted colloidal 
particles. We can control the interaction between polymer-grafted particles by adjusting the length, 
species, and the grafting density of the grafted polymers and the quality of the solvent [4-7]. Thanks 
to such usefulness, polymer-grafted colloidal particles can be adapted to a variety of applications and 
play an important role in industry. We study the effect of such grafted polymers on colloidal particles 
by Monte Carlo simulations using a simple model system. In the present article, we show that parti- 
cles interacting via spherically symmetrical repulsive potential, without any attraction, self-assemble 
and finally align in strings, and that this string-like assembly possesses a similarity to the percolation 
transition and critical phenomena. 

A pair of polymer-grafted colloidal particles interact not only with their rigid cores but also with 
their polymers outside the rigid cores. We can describe the interaction between the polymer-grafted 
colloidal particles by adding interaction potential of grafted polymers to hard core potential. For 
example, if the solvent is water, we graft block copolymers composed of a long hydrophobic block 
and a short hydrophilic block onto the colloidal particle onto which the end of the hydrophobic block 
is grafted. These grafted block copolymers form a melt brush covered with a thin swollen shell. Here a 
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Fig. 1 - The results of the SCF calculations on the interaction potential between a pair of the polymer-grafted 
colloidal particles with diameter <Ti for various polymerization index of grafted chains, N SG . We set the size of 
a segment 0.1<xi, the polymerization index of the matrix chains N S m = 5.0, the block ratio of the A-subchain of 
the grafted chains 0.998, the total volume of the segments of the grafted chains on both the colloidal particles 
\A6ttcti 3 , xab = 1-0, Xam = 1.0, Xbm = 0.0 where the indices A, B, and M indicate A segment, B segment, 
and the segment of the matrix chain respectively. At N S g = 20, these parameters correspond, for example, to 
the system obtained by grafting polypropylene with molecular weight « 10 5 onto colloidal particles 100 [nm] in 
diameter with the graft density « 1 [(nm) -2 ]. 



melt brush is defined as a brush without the solvent while a swollen brush is defined as a brush swollen 
by the solvent. 

Since two swollen brushes prevent themselves from overlapping with each other due to the ex- 
cluded volume effect, the interaction between a pair of swollen brushes is repulsive [8]. On the other 
hand, the interaction between two melt brushes is slightly attractive upon contact [8] since the chain 
conformation can relax when the brushes touch each other. The effect of the curvature of the grafting 
surfaces on the interaction between brushes can be evaluated using the so-called Derjaguin approxi- 
mation [9] as long as the radius of the curvature is small compared to the particle radius. However, if 
the radius of the curvature is the same order as the brush thickness, no good analytical result is known, 
and we have to rely on numerical self-consistent field (SCF) calculation [4-8]. 

First, we performed SCF calculations on two colloidal particles with diameter cr\ onto which A- 
B block copolymers with polymerization index N$g are grafted [10], where A and B subchains are 
hydrophobic and hydrophilic respectively. Figure [I] shows the calculated interaction potential as a 
function of r, the distance between the centers of the pair of the particles. The result indicates that the 
interaction potential approaches the square-step potential when Nsg increases while the total volume 
of segments of the grafted polymers is fixed. 

We approximate this potential by spherically symmetrical step potential with a rigid core, i.e. 
0(r) = oo (r < cr\), cp(r) = 6o ((T\ < r < cr 2 ), and cp(r) = (cr 2 < r), where cr 2 is the diameter of 
the outer core formed by the polymer brush and eo > 0. Note that the height of the potential step €q is 
dependent on the temperature because it originates from the conformational entropy. 

Owing to the step in the potential, the phase behavior of the system depends on temperature, which 
differs from the behavior of a rigid particle system [2, 3]. The probability that a pair of particles get on 
the step in the potential is given by the Boltzmann factor exp(-6o/^r). The phase behavior at extreme 
temperatures is easily understood though such a temperature region is out of the validity range of our 
model potential. 

At extremely high temperature the phase diagram of our system coincides with that of a system 
of rigid spherical particles with diameter <j\ since the step in the interaction potential is neglegible 
compared to the kinetic energy. 

At extremely low temperature, on the other hand, the step in the interaction potential becomes 
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extremely high compared with the thermal energy. Owing to the potential energy barrier the phase 
diagram also coincides with that of a system of rigid spherical particles with diameter cr 2 when the 
system volume, v, exceeds the close-packed volume of the outer core, vo. Then, the pressure jumps 
suddenly when the system is compressed across vo, which leads to a phase transition. 

In previous studies, a step potential with a rigid core (or one with both a rigid core and a square- 
well) was used as a model potential for cesium, cerium, water and so on to reproduce qualitatively the 
anomalies on their phase diagrams such as the anomalies in the behavior of the melting curves and the 
critical behavior of the phase diagram [11, 12]. In two dimensions the model system forms dimers, 
lamellae, and other structures [13-17]. In the present article we perform a series of simulations in 
order to obtain the complete phase diagram. 

In two dimensions Monte Carlo simulations based on the canonical ensemble have been carried out 
using the Metropolis algorithm [18]. Mersenne Twister is adopted as a random number generator [19]. 
In the initial state, N particles are arranged on a triangular lattice in a box of volume (area) v with 
periodic boundary conditions. At every simulation step a particle is picked at random and given a 
uniform random trial displacement within a radius of 0.2<ti. A Monte Carlo step (MCS) is defined 
as every N trial moves, during which each particle experiences one trial move on average. From now 
on dimensionless width of the potential step o- 2 /cri is always fixed at 2.0. From beginning to end of 
one simulation, the number of particles N, the dimensionless volume (area) v/vo (vo = No-2 2 V3/2 in 
2 dimensions), and the dimensionless temperature ksT /eq are fixed. After 5.0 x 10 6 MCS, by which 
the system relaxes to the equilibrium state except at low temperature and high density, we acquire data 
every 10 4 MCS and get 100 independent samples of particle configurations. In the present article, 
simulation results at N = 1200 are mainly shown. We have confirmed, however, that the physical 
properties of the simulation system are essentially unchanged even when we simulate a larger system 
with N = 4800. 

As a particular example, we show a snapshot of the system at N = 9900, v/vo = 0.7, ksT /cq = 0.1 
and 5.0 x 10 6 MCS in Fig.|3 We observe that the colloidal particles start to self-assemble and fi- 
nally align in strings. We call this effect "frogspawn effect" on account of the characteristic con- 
figuration [13-15]. With use of the molecular dynamics simulation based on the micro-canonical 
ensemble at N = 140, v/vo = 0.7, and internal energy per particle = 0.966o, which corresponds to 
ksT/eo « 0.036, we checked that this string-like assembly is not an artifact of the Monte Carlo simu- 
lation and actually observed the same string-like assembly. This configuration in 2 dimensions is also 
observed in a model system with cr 2 / cr\ = 2.0 [16] or 2.5 [17], where the system is cooled off from 
high temperature to low temperature, or in a model system that is composed of particles interacting 
via continuous repulsive potential [20] . 

To discuss the structure of the aggregation, including string-like one, we define string length and 
cluster size. First we define a "bond" as an overlap between a pair of particles. When bonds form a 
single network composed of M particles, we define this network as a "cluster with cluster size M". A 
particle that has more than two bonds in a cluster is defined as a "junction". 

In Fig. |21 we observe that clusters have different configurations of the network even if they have 
the same size. In order to quantitatively specify the network configurations of clusters, we introduce 
"string length". A string is defined as a sequence of particles, each of which has exactly two bonds. 
String length is defined as the number of particles along the string, including both ends, when the 
ends are junctions or particles with only one bond. We exclude a circular cluster composed of three 
particles from the definition of the strings though all the particles in such a cluster have two bonds. 

In order to check the dependence of the final particle configurations on the initial configurations, 
simulations with the same set of parameters have been performed from 4 different homogeneous and/or 
inhomogeneous initial configurations, which are perfect string-like assembly aligned in the x-direction, 
the close-packed configuration of the inner cores, etc. We calculate the string length distributions 
and the cluster size distributions for these 4 cases. These two distributions characterize the strings 
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Fig. 2 - The snapshot at v/vo = 0.7 and k B T/eo = 0.1. 

and clusters. The calculated distributions show that all the simulations starting from different initial 
configurations end in the same string-like structure, whose cluster size shows exponential distribution. 
This result indicates that the string-like configuration is the stable configuration. The exponential 
distribution of the cluster size indicates that the clusters are results of independent aggregation events 
of particles. This is just like the formation of clusters in percolation system in its non-critical region. 

In view of the internal energy we can understand the reason for the string-like assembly. At 
extremely low temperature the particles tend to get down the step in the potential to decrease internal 
energy. At v/vo < 1.0, however, the system does not have enough free volume to allow all the particles 
to get down the step. Consequently, at extremely low temperature and v/vo < 1.0, the particles form 
long strings since a string-like linear cluster takes lower internal energy than a concentrated bunching 
cluster does. As the temperature increases, the strings become short since the particles can climb the 
potential step more easily. Note that the particles do not form long strings at v/vo « 1.0 since the 
particles seldom overlap. At high density the particles compose only a few strings since the small free 
space allows very few particles to get down the step, and thus most particles are junctions. 

In order to determine the region of string-like configurations, we next calculate the number fraction 
of the particles that belong to any strings and the average string length. The average string length is 
shown in Figs. 0a) and (b). we observe that the average string length has a sharp peak at 0.1 < 
kBT/eo < 0.2 and 0.6 < v/vo < 0.7. On the other hand, the calculated number fraction (data not 
shown in this article) shows wide regions of string-like assembly, located at 0.0 < kBT/eo < 1.5 and 
0.3 < v/vo < 1.5. These results indicate that there are many short strings in these wide regions outside 
this peak. At high density regions, however, there are only a small number of strings since very few 
particles can get down the step. At v/vo > 1.0 and low temperature regions particles can seldom 
overlap and most particles keep isolated. Figure |5Jb) also indicates that the maximum average string 
length exceeds 100 at this peak, which suggests a divergence of the string length in that region. 

At N = 1200 and kBT/eo = 0.1 cluster size distribution (data not shown in this article) suddenly 
changes between v/vo = 0.6 and v/vo = 0.7. Atv/vo = 0.7 and 0.8 the distribution has only one 
peak at cluster size 2, whereas at v/vq = 0.6 and 0.5 a second peak appears around cluster size 1000. 
At v/vo > 0.7 a large majority of particles compose many small clusters in the system, whereas at 
v/vo < 0.6 most particles form one large cluster. In order to understand the sudden change in the cluster 
size distribution, we evaluate the probability of the occurrence of a large cluster that bridges both sides 
of the system, i.e. a percolated cluster [21]. In Fig.@Ja), we show the occurence probability of the 
percolated cluster that is obtained from the 100 independent samples at each point of the (k B T/eo)- 
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Fig. 3 - The average string length at N = 1200. Blue regions represent short strings and red ones long strings. 
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Fig. 4 - Occurrence probability of percolated clusters, (a), and its details around the boundary, (b), for the same 
systems as Fig.|3] 



(v/vo) diagram. We observe that the percolated clusters suddenly appear when the system crosses 
a sharp boundary around v/vo ~ 0.7-1.0. To determine the location of this boundary, we show in 
Fig. |4jb) a similar picture to Fig. |Ua), but with a finer resolution. We can see that the occurrence 
probability of the percolated clusters changes abruptly on a sharp boundary, which suggests a similarity 
of this phenomenon to the usual percolation phenomena [21]. 

The boundary in Fig.|4{a) asymptotically approaches to the line v/vo = 1.0 at high temperature, 
since at v/vq > 1.0 the system has enough free volume to allow all particles to get down the step and 
the particles do not frequently overlap. On the other hand, at low temperature percolation transition 
occurs at v/vo ~ 0.65, between Alder transition of outer cores and Alder transition of inner cores. 

Now we study the relationship between our string-like assembly and the percolation phenomena. 
We calculate the cluster size distribution along the boundary curve in Figs. @Ia) and (b) (data not 
shown in this paper). The result shows that the cluster size distribution is a power law function of 
the cluster size with the exponent about 1.89. This exponent 1.89 corresponds to the value of the 
size distribution of percolated clusters on a 2-d triangular lattice, which equals 2.0 [21]. Out of the 
percolation transition region, the cluster size distribution becomes an exponential function, as was 
already mentioned, which is also similar to the usual percolation phenomena. 

Figures 0a) and (b) also show an anomaly, i.e. a divergence of the string length, around the 
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Fig. 5 - Phase diagrams in 2 dimensions at N = 1200, (a), and in 3 dimensions at N = 4000, (b). •: Solid phase. 
□ : Fluid phase or coexistence of both phases. Dotted lines correspond to the Alder transition. 




Fig. 6 - The average string length, (a), and occurrence probability of percolated clusters, (b), in 3 dimensions. 
N = 4000. 

boundary between the percolated region and the non-percolated region. This relation between string- 
like configuration and percolation transition indicates a similarity between our string-like assembly 
and the critical phenomena. Due to the poor resolution of the data, however, we cannot evaluate the 
critical exponent. 

By calculating the mean square displacement (MSD) of the particles during 10 6 MCS, we deter- 
mine regions of a solid phase and draw a phase diagram (see Fig.|5Ja)). We regard the system as a solid 
if the MSD value is less than its value at ksT /eq = 0.01 and v/vo = 1.3. This is because Alder transition 
occurs at v/vo ~ 1.3 in 2D hard particle systems [2, 3]. At high temperature regions Alder transition 
of inner cores occurs at 0.3 < v/vo < 0.4, whereas Alder transition of outer cores, which is expected 
at v/vo ~ 1.3, does not. At low temperature regions only Alder transition of outer cores emerges. At 
intermediate temperature regions the solid melts when not only the volume but also the temperature 
increases. These results correspond to the phase behavior discussed in the introduction. The region 
where our string-like assembly like Fig. Eland percolation transition occur, i.e. 0.1 < k B T/eo < 0.2 
and 0.6 < v/vq < 0.7, is located in a solid phase between two Alder transitions. We can see that our 
string-like assembly like Fig.|2|is an amorphous solid. 

In three dimensions Monte Carlo simulations have also been performed according to the same 
simulation condition as 2D simulations, except that the particles are arranged on fee lattice in the 
initial state and we acquire data after 10 6 MCS. Figures EJb),[6ja), and|6{b) show a phase diagram, the 
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average string length, and occurrence probability of percolated clusters in 3 dimensions respectively. 
At the percolation transition the cluster size distribution becomes a power law function of the cluster 
size with the exponent about 2.1, which is different from the value 1.89 for 2D case. These results 
qualitatively agree with the ones in 2 dimensions, which indicates that our discussion on 2D systems 
also applies to 3D systems. Our results in both 2 and 3 dimensions show that the average string 
length diverges around the region where the melting transition line and the percolation transition line 
cross. A similar behavior is found in Ising spin system where the percolation transition line and the 
order-disorder line meet at the critical point [21]. 

In conclusion, using Monte Carlo simulations we studied phase behavior of 2 and 3 dimensional 
systems of hard particles with step repulsive interacting core. The simulation results show string- 
like assembly of particles between the two Alder transition densities, each of which corresponds to 
the inner hard core and the outer hard core. We found that this string-like assembly is related to the 
percolation phenomena. We can also observe the same string-like assembly like Fig.|2|in 2 dimensional 
systems at cr 2 /cri = 1.9, 2.1 and 2.5 though not at o- 2 /cri = 1.1, 1.5, 3.0, 3.5 and 10.0, which will be 
discussed in our forthcoming paper [10]. 
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